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Abstract 

Background: Historically, the Chittagong Hill Tracts (CH~0 of Bangladesh was considered hyperendemic for malaria. 
To better understand the contemporary malaria epidemiology and to develop new and innovative control 
strategies, comprehensive epidemiologic studies are ongoing in two endemic unions of Bandarban district of CHT. 
Within these studies entomological surveillance has been undertaken to study the role of the existing anopheline 
species involved in the malaria transmission cycle throughout the year. 

Methods: CDC miniature light traps were deployed to collect anopheline mosquitoes from the sleeping room of 
the selected houses each month in a single union (Kuhalong). Molecular identification was carried out for available 
Anopheles species complexes. Circumsporozoite proteins (CSP) for Plasmodium falciparum, Plasmodium wVax-210 
(Pv-210) and Plasmodium wVax-247(Pv-247) were detected by Enzyme-linked immunosorbent assay (ELISA) from the 
female anopheline mosquitoes. To confirm CSP-ELISA results, polymerase chain reaction (PCR) was also performed. 

Results: A total of 2,837 anopheline mosquitoes, of which 2,576 were female, belonging to 20 species were 
collected from July 2009 -June 2010. Anopheles jeypohensis was the most abundant species (18.9%), followed by An. 
vagus (16.8%) and An. kochi (14.4%). ELISA was performed on 2,467 female mosquitoes of 19 species. 15 (0.6%) 
female anophelines belonging to eight species were found to be positive for Plasmodium infection by CSP-ELISA. 
Of those, 1 1 (0.4%) mosquitoes were positive for P. falciparum and four (0.2%) for Pv-210. No mosquito was found 
positive for Pv-247. An. maculatus (2.1%, 2/97) had the highest infection rate followed by An. umbrosus (1.7%, 2/1 15) 
and An. barbirostris (1.1%, 2/186). Other infected species were An. nigerrimus, An. nivipes, An. jeypohensis, An. kochi, 
and An. vagus. Out of 1 1 P. falciparum CSP positive samples, seven turned out to be positive by PCR. None of the 
samples positive for Pv-210 was positive by PCR. In terms of abundance and incrimination, the results suggest that 
An. maculatus, An. jeypohensis and An. nivipes play important roles in malaria transmission in Kuhalong. 

Conclusion: The findings of this study suggest that even in the presence of an insecticide impregnated bed-net 
intervention, a number of Anopheles species still play a role in the transmission of malaria. Further investigations are 
required to reveal the detailed biology and insecticide resistance patterns of the vector mosquito species in 
endemic areas in Bangladesh in order to assist with the planning and implementation of improved malaria control 
strategies. 
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Background 

In Bangladesh, malaria is most endemic in the highland 
areas bordering India and Myanmar. Of the 13 malaria 
endemic districts in Bangladesh, the Chittagong Hill 
Tracts (CHT) are composed of three districts; Ranga- 
mati, Bandarban and Khagrachhari, which contribute 
80% of the country's total malaria case burden [1]. Plas- 
modium falciparum and P. vivax are the two main para- 
sites in CHT [1,2] although some sporadic reports of P. 
malariae and P. ovale have been reported recently [3,4]. 
Malaria incidence is seasonal in Bangladesh where the 
warm and wet months of May- October define the peak 
season and the dry and cooler months of November- 
April define the off season [5]. 

To date, 35 anopheline species have been reported in 
Bangladesh [6]. Of those, four species, An. baimaii {dims 
D), An. minimus, An. sundaicus and An. philippinensis 
were incriminated as malaria vectors during the Malaria 
Eradication Programme (MEP) of the 1960s. In the early 
1990s three additional species {An. aconitus, An. annu- 
laris and An. vagus) were incriminated following three re- 
spective outbreaks in flood plain areas of the country 
[7,8]. However, due to loss of forest habitat, the density of 
two primary vectors' {An. baimaii and An. minimus) were 
dramatically reduced in several sporadic entomological 
investigations carried out by the Malaria and Parasitic Dis- 
ease Control Unit (M&PDC) (personal communication 
with N. Chaudhury, Senior Entomologist, M&PDC). 

An entomological investigation carried out in three 
malaria endemic border areas of the country during the 
peak of the 2009 transmission season found seven 
Anopheles species harboring Plasmodium based on 
ELISA. Except for An. philippinensis and An. vagus five 
other species {An. karwari, An. maculatus, An. barbiros- 
tris, An. nigerrimus, and An. subpictus), not previously 
incriminated in Bangladesh, were found to be infected 
[8]. From that study it was concluded that in the absence 
of recognized primary malaria vectors', other Anopheles 
species might play a significant role in continued malaria 
transmission in Bangladesh. 

In collaboration with the Johns Hopkins Malaria Re- 
search Institute (JHMRI), icddr,b initiated a study to 
understand the epidemiology of malaria transmission in 
two unions (Kuhalong and Rajbila) of Bandarban district 
in mid 2009. The method has been described elsewhere 
[9]. The distribution of the anopheline fauna and the 
species' roles in the malaria transmission cycle in the 
study areas during both the wet and dry season remain 
the main focus of the current entomological investiga- 
tions. The aim of this paper is to report the faunal diver- 
sity and Plasmodium infection status of anopheline 
species that have been collected in the first year (July 
2009- June 2010) of the entomological surveillance in 
Kuhalong. 



Methods 

Study area 

Kuhalong (22° 12' 45" N, 92° 9' 35" E) is located adjacent 
to the Bandarban town. Kuhalong was divided into 12 
clusters, which had roughly equal numbers of house- 
holds as part of the epidemiologic study. A description 
of Kuhalong and its clusters have been provided earlier 
[9]. In brief, the union has a total area of 79 sq km with 
a population of more than 11,390. Kuhalong is a hilly 
and forested area with an average elevation of 80 meters 
ranging up to 152 meters. Some other important fea- 
tures of Kuhalong include the presence of intricate 
branches of rivers, several streams, marshy lands and 
plant monoculture (teak and rubber). The land for plant 
monoculture was created by destroying natural forest. 
99.5% of the inhabitants of the area have bed net cover- 
age (untreated/insecticide treated). But when asked, 10% 
of the respondents admitted that they had not used a 
bed net the previous night (unpublished data). 

Collection of Anopheles mosquitoes 

CDC miniature light trap (model 512, John W. Hock 
Inc, USA) was deployed for indoor mosquito collection 
from the sleeping room. Traps were placed for 12 hours 
(6 pm to 6 am). During the wet season 100 houses were 
selected randomly, based on elevation and were trapped. 
This was reduced to 50 houses in the off (dry) season. 
Later, when trapping started in the second Union (Raj- 
bila), the number was fixed to 60 houses at each union 
(5 houses in each of 12 clusters) due to the available 
resources and personnel. Each house was later trapped 
once in a month throughout the season in order to ob- 
tain longitudinal information. 

Mosquito sample preparation 

Mosquitoes were collected from the light traps using a 
battery-powered aspirator and were subsequently killed 
by chloroform. Anopheline mosquitoes were sorted and 
preserved temporarily in plastic tubes. After identifying 
the species at the field office in Bandarban each mos- 
quito was preserved in a separate plastic vial (1.5 ml) la- 
belled and capped with silica gel and cotton. 

Circumsporozoite protein (CSP) ELISA 

The ELISA methods were described previously else- 
where [10,11]. ELISA distinguished between circumspor- 
ozoite protein (CSP) of P. falciparum and two distinct 
polymorphs of P. vivax: Pv-210 and Pv-247. In each test, 
field caught male Anopheles spp. was used as the nega- 
tive control. The positive controls and monoclonal anti- 
body (MAB) for the ELISAs were obtained from the 
Centers for Disease Control and Prevention, Atlanta, 
USA. The optical density (OD) was measured at 410 nm 
in a Bio-Rad ELISA plate reader, 60 minutes after adding 
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the substrate [8]. A cut-off value at least twice the mean 
OD of the negative controls was considered as positive. 
ELISA positive samples in a given day were repeated the 
next day to confirm the result. 

PCR 

Remaining mosquito lysates were used for DNA extrac- 
tion using the QiaAmp DNA Minikit (Qiagen, Valencia, 
CA), following the tissue extraction protocol. Nested 
PCR was performed to detect P. falciparum and P. vivax 
infection from extracted mosquito DNA using the pri- 
mers described by Snounou et al. [12]. PCR conditions 
were modified to adjust for low concentration of DNA 
from mosquito samples. 2 \iL of the amplicon from the 
first step was used as template in the second step. 

To confirm taxonomic identification of some Anoph- 
eles species, a confirmatory PCR was made available, 
based on availability of the established protocols. Proto- 
cols used were Phuc et al. [13] for species belonging to 
An. minimus and Myzomia Series, two different 



protocols of Walton et al. [14,15] for Annularis and 
Maculatus groups, Huong et al. [16] for An. dims sibling 
species, and Goswami et al. [17] for the An. culicifacies 
complex, respectively. Protocols were followed with 
some modifications. DNA amplifications were done on 
S1000® Thermal Cycler (Bio-Rad Laboratories, Inc., 
Hercules, CA) and visualized under UV illumination 
after electrophoresis on ethidium bromide-stained 1.5% 
agarose gels along with invitrogen® 100 base pair (bp) 
molecular mass marker (Life technologies, NY, USA). 
Representative samples (DNA from homogenates) were 
sent for sequencing to the JHMRI. The list of primers is 
given in Additional file 1. 

Sequencing of mosquito specimens 

The variable internal transcribed spacer 2 region (ITS2) 
was amplified with primers from the flanking 5.8 S and 
28 S genes; ITS2A (5-TGT GAA CTG CAG GAC 
ACA T-3 ) and ITS2B (5-ACC CCC TGA ATT TAA 
GCA TA-3') [18]. Because sequencing reactions with 






Figure 1 A satellite image over study area showing location of 268 trapping houses. Here, positive means houses that had CSP-positive 
mosquitoes and negative means houses with CSP-negative mosquitoes or without female anopheline mosquitoes. 
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ITS2A frequently failed, some samples were amplified 
and sequenced using the novel primer ITS2B1 (5'-GTC 
CCT ACG TGC TGA GCT TC-3'). This primer binds 
further downstream in the 28 S gene; such that the 3' 
portions of products immediately upstream of the 
ITS2B binding site, could be sequenced. Each 25 uL re- 
action contained IX PCR buffer, 200 uM each dNTP, 30 
pmol each primer, 2 units Taq polymerase, and 1.0 uL 
DNA template. Products were amplified in a thermocy- 
cler (MJ Research, Watertown, MA, USA) using the fol- 
lowing conditions: 2 minute initial denaturation at 94°C, 
40 cycles of 30 seconds at 94°C, 30 seconds at 50°C, 
and 40 seconds at 72°C, and a 10 minute final extension 
at 72°C. PCR product size ranged from 480 bp to 
847 bp. 5 uL of each PCR product was subjected to 
electrophoresis on a 2% agarose gel, stained with eth- 
idium bromide, and visualized under UV illumination. 
The remainder of each successful PCR reaction was 
purified using a Qiaquick PCR prep kit (Qiagen, Valen- 
cia, CA). Products were sequenced in both directions 
using dye terminator chemistry on a 3730x1 DNA 
Analyzer (Applied Biosystems, Foster City, CA) at the 
Johns Hopkins University School of Medicine Sequen- 
cing and Synthesis Facility. 



Results 

Collections were made during 641 trap nights from July 
2009 to June 2010 in households throughout Kuhalong 
(Figure 1). A total of 2,837 Anopheles mosquitoes (2,576 
female and 261 male) were caught (4.4 mosquitoes/trap 
night; std. error 0.21). Among the collected female ano- 
phelines, 20 species were confirmed based on taxonomic 
characteristics, adjusted by PCR for available species com- 
plexes. A representative half of the specimens belonging 
to An. philippinensis based on morphological characteris- 
tics were later confirmed as An. nivipes by PCR followed 
by sequencing. Similarly the identity of all An. pallidus 
samples (approximately 30 specimens) was revised to An. 
nivipes after molecular confirmation. Therefore, all indi- 
vidual specimens originally morphologically identified as 
An. philippinensis and An. pallidus were listed as An. 
nivipes in our final data set. A number of An. minimus 
were molecularly identified as An. varuna, whereas the 
remaining An. minimus specimens were confirmed as An. 
minimus subspecies A. Three subspecies of An. culicifacies 
(B, C and E) were confirmed, with dominance of subspe- 
cies B. A strong concordance was observed between mor- 
phological and moleculer detection for An. baimaii and 
An. maculatus group. 





Figure 2 Overall contribution of female Anopheles species in Kuhalong. 
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Three species were nearly co-dominant. An. jeyporien- 
sis (18.9%) was the most commonly captured species, 
followed by An. vagus (16.8%) and An. kochi (14.4%) 
(Figure 2). The relative abundance and dominance of 
anopheline species varied temporally throughout the 
year (Figure 3). Slightly more female anophelines were 
caught in the non-malaria dry season (n = 1391, mean 
4.3 mosquitoes/trap night) than peak wet transmission 
season (n = 1185, mean 3.7 mosquitoes/trap night). 
However, this seasonal difference in mosquito collection 
rates was not statistically significant (p>0.05). The high- 
est diversity of species was found in July 2009 (16) and 
lowest in September 2009 and April 2010 (11). An. 
nivipes appeared to be the dominant species during peak 
malaria season. However, in the off season, An. kochi 
was the dominant species. In both seasons, An. jeypor- 
iensis and An. vagus were the second and third most nu- 
merous species, respectively. The density of An. nivipes 
complex was very low in the off season. An. maculatus, 
An. nigerrimus, An. umbrosus and An. barbirostris 
appeared to retain a moderate density throughout the 
year (Figure 4). 

CSP-ELISA was performed on 2,467 female anopheline 
mosquitoes identified to 19 species. The untested mos- 
quitoes (n = 109), including a single specimen of An. tur- 
khudi, were kept as voucher specimens. Fifteen 
mosquitoes collected from 13 houses and belonging to 
eight species were Plasmodium-positive by CSP-ELISA. 
Thus, the overall infection rate was 0.6% (15/2467). By 
species, the highest infection rate was observed in An. 
maculatus (2.1%) followed by An. umbrosus (1.7%), An. 
barbirostris (1.1%), An. nigerrimus (1%), An. nivipes 
(0.8%), An. jeyporiensis (0.6%), An. kochi (0.5%), and An. 
vagus (0.2%), respectively (Table 1). 



Among the CSP-positive anophelines, 11 (0.4%) mos- 
quitoes belonging to six species were found positive for 
P. falciparum and four (0.2%) mosquitoes belonging to 
three species were found positive for Pv-210 (Table 2). 
None of the mosquitoes were positive for Pv-247 or had 
mixed Plasmodium infection. P. falciparum positive 
Anopheles species included An. barbirostris, An. jeypor- 
iensis, An. kochi, An. maculatus, An. nigerrimus, and An. 
nivipes. Pv-210 positive species included An. nivipes, An. 
umbrosus and An. vagus. 

The highest monthly infection rates for Plasmodium 
were observed in October 2009 (2%) and November 
2009 (1.3%). Plasmodium-infected mosquitoes were not 
detected in September or December of 2009 and April- 
June 2010 (Figure 5). In the dry season three sporozoite 
positive mosquitoes were caught in November 2009 and 
two in March 2010. Although the infection rate in peak 
season (0.8%) was twice that of off season (0.4%), this 
seasonal difference in mosquito infection rates was not 
statistically significant (p>0.05). 

Of the 15 Plasmodium-positive mosquitoes, seven of 
these were positive for the presence of P. falciparum by 
nested PCR out of 11 that tested positive by CSP-ELISA. 
However, none of the Pv-210 CSP positive samples were 
positive by PCR. Species that had both CSP and PCR 
positive samples included five species: An. maculatus, 
An. jeyporiensis, An. nivipes, An. barbirostris and An. 
nigerrimus (Table 2). 

Discussion 

The Anopheles species diversity is very high in Kuhalong. 
This diversity is even higher than what had been 
observed in a previous study in Matiranga, a sub-district 
of CHT, where 15 anopheline species were reported in a 
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Figure 3 Monthly frequency of five prevalent female Anopheles species in Kuhalong. 
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Figure 4 Contribution of eight prevalent female Anopheles species according to transmission season. 



peak transmission season, of which seven were positive 
for CSP [8]. The greater number of anopheline species 
collected in Kuhalong might be due to the cumulative 
year-long trapping effort, as several species were most 
frequently sampled during the low transmission season 
(Figures 3 and 4). This highlights the importance of off 
season trapping. 

Of the eight Plasmodium-infected anopheline species 
reported in this study, An. jeyporiensis and An. kochi 
have never been previously incriminated as vectors in 
Bangladesh. However, An. jeyporiensis and An. kochi 
have been considered to be potential vectors in the foot- 
hills of Assam, a north-eastern Indian state [19]. There 
is also a report of sporozoites in An. jeyporiensis made 
in 1944 by Macan from mosquitoes collected on the 
Myanmar-Bangladesh border [20]. The high relative 
abundance of An. jeyporiensis and its infection con- 
firmed by CSP followed by PCR indicates that this spe- 
cies likely plays an important role in the malaria 
transmission in the area. 

Although An. maculatus is considered an important 
vector in Peninsular Malaysia and Southern Thailand 
[21,22], it was thought to be a less important species in 
Bangladesh earlier. This species was first found to be 
CSP infected in Bangladesh in recent times on the basis 
of a single positive mosquito out of seven tested [8]. Our 
results re-emphasize its potential importance as a vector 
in this region. 



Despite expectations based on prior observations, no 
An. baimaii or An. minimus s.l. were found CSP-positive 
during this study in Kuhalong. Further study is required, 
but these observations may indicate that these species 
are less abundant, are locally less competent or that this 
reflects a foraging effect/trapping bias. Similarly, no An. 
culicifacies, a major malaria vector in some parts of 
India and Sri Lanka [23,24], were found to be infected 
with Plasmodium in Kuhalong. It is worth noting that 
An. culicifacies is primarily represented by species B in 
Kuhalong, a species considered to be a non-vector or a 
poor vector [24,25], and it is common in the eastern 
states of India sharing borders with Bangladesh [26]. 
This supports the idea that this species does not have a 
vectorial role in Kuhalong. Apart from An. vagus, two 
other potential secondary malaria vectors in Bangladesh, 
An. aconitus and An. annularis [27,28], were rarely col- 
lected in Kuhalong and all were negative for Plasmo- 
dium infection by CSP ELISA. 

An. philippinensis has historically been considered the 
main vector of malaria in the vast plain areas of Bangla- 
desh, but has been considered rare in the CHT [7]. The 
findings of the present study, in accordance with other 
recent observations [8], initially suggested that this spe- 
cies has adapted to the high land areas and has 
expanded its role in the transmission of malaria. How- 
ever, An. nivipes is morphologically very similar with An. 
philippinensis and often confused with the latter. The 
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Table 1 Prevalence of female anopheline mosquitoes and 
their CSP positive rates 
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adults of these two species are differentiated by a single 
characteristic related to a pre-sector dark mark on the 
wing vein [29] . In the Assam State of India, recent stud- 
ies document the role of the An. philippinensis-nivipes 
complex in malaria transmission [19,30] and later con- 
firmed the presence of P. falciparum DNA in two An. 
nivipes specimens, that were molecularly identified [31]. 
The presence of P. falciparum specific CSP by ELISA 
and later parasite DNA in An. nivipes in the present 
study support the previous observations in India. These 
two species have never been differentiated in Bangladesh 
until this report. 

Representatives of An. barbirostris and An. nigerrimus 
species were found to be positive for the presence of 
CSP, which was also seen in previous observations made 
in the CHT [8]. Previously in this region, An. barbiros- 
tris was confirmed as a vector in Sri Lanka based on 
CSP-ELISA [32]. 

Depletion of the forest habitats, preferred by An. bai- 
maii and An. minimus and conversion to agricultural 
use may have allowed members of the An. nivipes to ex- 
tend their range into highland areas. It is unknown how 
these species interact, whether it is the habitat that 
drives the relative abundance of these species or whether 
there is active competition and replacement as well. 



Thus, the bionomics of anopheline species requires fur- 
ther study and will reveal valuable information necessary 
for future control strategies. 

As a method, CSP-ELISA has several advantages al- 
though it is thought to be less sensitive, particularly 
when low numbers of sporozoites are present [33]. 
Hence, a false positive result is also possible due to the 
presence of CSP in body parts other than salivary glands 
[33,34] and also due to the presence of bovine or swine 
blood [35]. On the other hand, PCR is more sensitive 
than ELISA and can detect fewer parasites [36], but is 
not sporozoite specific and is recommended to recon- 
firm CSP-ELISA results [37]. None of the Pv-210 CSP 
positive samples in the study tested positive by PCR. A 
low concentration of DNA obtained from CSP homoge- 
nates could be another reason for negative PCR results 
along with CSP false positive samples. 

In terms of incrimination by overall abundance and in- 
fection rates, An. maculatus, An. jeyporiensis and An. 
nivipes all are likely to have important roles in malaria 
transmission in Kuhalong. This is not to say that other 
species do not also contribute significantly to malaria 
transmission in Kuhalong. 

The National Malaria Control Programme (NMCP) of 
Bangladesh initiated a programme in 2007 with financial 
support from the Global Fund. This programme includes 
diagnosis and treatment, distribution of long-lasting in- 
secticidal bed net (LLIN) and re-treatment of 
insecticide-treated nets (ITNs). As of 2010, the NMCP 
has treated/retreated more than 2 million bed nets with 



Table 2 Summary information of CSP positive Anopheles 
female mosquitoes from Kuhalong, Bandarban 
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Pf 


Pf 


0 


09-0459 


An. umbrosus 


Pv-210 


0 


1 


09-0775 


An. jeyporiensis 


Pf 


0 


0 


09-0785 


An. umbrosus 


Pv-210 


0 


1 


09-0826 


An. barbirostris 


Pf 


0 


0 


10-1147 


An. jeyporiensis 


Pf 


Pf 


0 


10-1700 


An. kochi 


Pf 


0 


0 


10-1917 


An. nivipes 


Pv-210 


0 


1 



a Pf=P. falciparum, Pv-210 = P. v/Vax-210. 
b Pf = P. falciparum{+ve), 0 = negative. 
c 0 = unfed, 1 = blood fed. 
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insecticides and distributed almost 1.7 million LLINs 
throughout the malaria endemic areas of the country 
[38]. However, no prominent work has been reported so 
far to monitor the resistance patterns of the anopheline 
mosquitoes against the insecticides used in these control 
efforts. The high diversity of anophelines, may increase 
the risk of insecticide resistance and could alter biting 
behaviours which may circumvent some of the protec- 
tion from bednets. This high species diversity is typical 
of Southeast Asia, in contrast to many parts of Africa 
where only a few species are major contributors to Plas- 
modium transmission. 

Conclusions 

The findings of this study illustrate that even in areas of 
Bangladesh where ITNs and LLINs have been deployed as 
malaria interventions, a number of diverse Anopheles spe- 
cies still play a role in the transmission of malaria. Further 
study of the bionomics, ecology, and insecticide resistance 
of these species is necessary to understand the transmis- 
sion biology and to provide the information required for 
the development of evidence -based control programmes 
in Bangladesh. This insight may be useful to other coun- 
tries similarly facing high burdens of malaria. 

Additional file 



Additional file 1: List of primers. 
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